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METHOD AND DEVICE FOR FITTING SURFACE TO POINT GROUP, 
MODELING DEVICE , AND COMPUTER PROGRAM 



This application is based on Japanese patent 
5 applications Nos. 2000-090629, 2000-197778 and 2000-197779 
filed in Japan, the contents of which are hereby 
incorporated by reference. 



BACKGROUND OF THE INVENTION 

10 1. Field of the Invention 

The present invention relates to a method and a 
device for fitting a surface to a point group, a modeling 
device and a computer program. The present invention is 
useful in the field of computer graphics such as 

15 generation of a three-dimensional model, etc. 
2. Description of the Prior Art 

In recent years, the three-dimensional CG (three- 
dimensional computer graphics) technique is often used in 
the fields of movie, game and so on. In the three- 

20 dimensional CG, a three-dimensional model, light source, 
etc. are located and moved in a virtual three-dimensional 
space; therefore, the three-dimensional CG is high in the 
degree of freedom of expression. 

Conventionally, a non-contact three-dimensional 

25 measurement device employing the light -section method has 
been put into practical use. Three-dimensional data of an 
object can be relatively easily obtained by using such 
non- contact three-dimensional measurement device. However, 
various problems have been detected in using the three - 

30 dimensional data obtained by means of the non-contact 



three-dimensional measurement device as it is for three- 
dimensional CG. For example, it is necessary to reduce 
the amount of data obtained by the measurement by a 
process of thinning the data, which is complicated and 
5 time-consuming. 

To solve the above problems , there has been proposed 
a method wherein a standard model is prepared for an 
object and the standard model is modified according to 
three-dimensional data of the object obtained by an actual 

10 measurement of the object (Japanese Unexamined Patent 
Publication No. 5-81377). 

In the above-mentioned method, three-dimensional 
form information of the three-dimensional data obtained by 
the measurement, i.e., a point group in the three - 

15 dimensional space, is used as an object for fitting, and a 
surface of the standard model is fitted on the three- 
dimensional point group. According to the method, it is 
possible to obtain three-dimensional data or a three- 
dimensional model free from defects even in the case where 

20 the three-dimensional data originally had defects. 

However, a satisfactory result is not yet achieved 
by the above method, since, in modifying the surface 
simply according to the point group obtained by the 
measurement, the surface is unavoidably fitted to some 

25 false points included in the point group, which are low in 
reliability. 

Specifically, as shown in Fig. 19 A, in the case 
where the points P3 to P6 among the six points PI to P6 
are low in reliability, the surface S is modified 
30 according to all the points PI to P6. The surface Sa 



obtained by the modification is not always satisfactory or 
proper . 

Further, as shown in Fig. 19B, in the case where 
points P3 to P6 , which are low in reliability, are ignored 
5 and reliable points PI and P2 are used, the modification 
itself cannot be fully achieved with failing to reflect 
and maintain a form of an object represented by the point 
group. The surface Sb obtained by the incomplete 
modification is hardly a proper representation of the form 
10 of the object. 

The three-dimensional model obtained by the above- 
mentioned method is frequently utilized in the field of 
animation. In the field of animation, not only a skin 
model which is external and visible for human eyes, but 
15 also a skeleton model or a muscle model (muscle 
information) which moves the skin model as lying 
thereunder is utilized. 

However, since the three-dimensional data obtained 
by measurement of an object corresponds to the skin model, 
20 it is impossible to generate the skeleton model or the 
like from the three-dimensional data as it is. 

Further, in the case of adopting modification 
information of a skin model as it is for generating a 
skeleton model based on the three-dimensional data of the 
25 object, it is possible that the skeleton model is 

improperly modified due to the difference between the skin 
model and the skeleton model . 

Although a general form of an object can be 
reproduced according to the above method, it is difficult 
30 to reproduce topical characteristics of a form of the 



object . 

For example, in the case of making a three- 
dimensional model of a human head according to the above 
method, although general form of the resulting three - 
5 dimensional model conforms to an object for measurement, 
i.e., the human head, parts that influences much for 
expression of the face such as eyes, mouth and the like do 
not satisfactorily conform to those of the object. Thus, 
there have been problems relating to the representation of 
10 delicate expressions of human. 

SUMMARY OF THE INVENTION 
An object of the present invention is to prevent 
improper modification or erroneous modification that is 
15 caused by a point of low-reliability and to perform proper 
fitting even when a point of low- reliability is included 
in a point group. 

Another object of the present invention is to 
generate a form model such as a skin model, skeleton model 
20 and the like by modifying a model for skin, a model for 

skeleton and the like with fitting using measurement data 
obtained by measuring an object. 

Another object of the present invention is to have 
parts of a three-dimensional model such as a corner of an 
25 eye or a corner of a mouth conformed to that of an object 
without topical improper modification. 

According to one embodiment of the present invention, 
there is provided a method for fitting a surface to a 
point group by using a computer which comprises a first 
30 step of judging reliability of the point group and a 



second step of changing the method for fitting the surface 
to the point group based on the judgment result obtained 
in the first step. 

Preferably, the first step may comprise judging 
5 reliability of each of points composing the point group, 
and the second step may comprise changing the method for 
fitting the surface to the point group by varying weight 
of each of the points based on reliability of each of the 
points . 

10 According to another embodiment of the present 

invention, there is provided a processor comprising an 
obtaining section for obtaining original data generated by 
measurements, a first modifying section for modifying, 
based on the original data, a first standard model which 

15 has been prepared separately from the obtaining of the 

original data and a second modifying section for modifying 
a second standard model based on the result of 
modification of the first standard model. 

Preferably, the first standard model may be a model 

20 for skin, the second standard model may be a model for 
skeleton or a model for muscle and the original data may 
be data whose object is a surface of a human head. 

According to a more preferable embodiment of the 
present invention, there is provided a device for 

25 modifying a surface based on a three-dimensional point 
group comprising a selecting section for selecting a 
plurality of partial areas from the point group and a 
modifying section for modifying the surface based on the 
point group of each of the selected partial areas . 

30 These and other characteristics and objects of the 
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present invention will hereinafter be described in detail 
by the following preferred embodiments. 



BRIEF DESCRIPTION OF THE DRAWING 
5 Fig. 1 is a block diagram showing the modeling 

device according to an embodiment of the present invention. 

Fig. 2 is a flow chart showing the general process 
of the modeling device. 

Fig. 3 is a flow chart showing the modification 
10 process. 

Fig. 4 shows an example of the standard model. 
Fig. 5 shows an example of the obtaining of three- 
dimensional data from an object. 

Fig . 6 shows an example of rough alignment . 
15 Fig. 7 shows an example of the process of extraction 

of outlines and characteristic points . 

Fig. 8 shows an example of a surface of a standard 
model and a point of three-dimensional data. 

Fig. 9 shows an example of a virtual spring for 
20 preventing a standard model from improper modification. 

Fig. 10 shows an example of the method of obtaining 
three-dimensional data and reliability data of an object. 
Fig. 11 shows an example of extraction of partial 

areas . 

25 Fig. 12 shows an example of the method of performing 

modification of a standard model using a partial area. 

Fig. 13 shows relationship between a model for skin 
and a model for skeleton with respect to construction 
points and control points. 

30 Fig. 14 is an example of a method of correcting a 
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control point. 

Fig. 15 is another example of the method of 
correcting a control point . 

Fig. 16 is an example of generation of a three- 
5 dimensional model by modifying a model for skin and a 
model for skeleton. 

Fig. 17 is a flow chart showing a process of 
generation of a three-dimensional model by modifying two 
models . 

10 Figs. 18A and 18B are examples of fittings of a 

surface S to a point P. 

Figs 19A and 19B are examples of conventional 
fittings of a surface S to a point P. 



15 DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig. 1 is a block diagram showing a modeling device 

1. 

The present embodiment illustrates an example of 
modifying a standard model based on measurement data 

20 (three-dimensional data or two-dimensional data) of a 
human head, i.e., an example of generating a three- 
dimensional model of a human head by fitting. 

The modeling device 1 generates a three-dimensional 
model such as a skin model and a skeleton model by 

25 modifying a standard model such as a model for skin, a 
model for skeleton and the like . Explanations for the 
generation of the skin model and the generation of the 
skeleton model are given below in this order. 
(Generation of Skin Model) 

30 Referring to Fig. 1, the modeling device 1 comprises 



a processor 10, a magnetic disk unit 11, a media drive 
unit 12, a display 13, a key board 14, a mouse 15, a 
three-dimensional measuring unit 16 and so on. 

The processor 10 comprises a CPU, a RAM, a ROM, a 
5 video RAM, an input and output port, various controllers 
and so on. The CPU executes programs stored in the RAM, 
the ROM, etc. , so that various functions are realized on 
the processor. 

Stored in the magnetic disk unit 11 are an OS 

10 (Operating System) , a modeling program PR for generating a 
three-dimensional model ML, other programs, a standard 
model (standard model data) DS, three-dimensional data 
(three-dimensional measurement data) DT, reliability data 
DR which indicate reliability of the three-dimensional 

15 data DT, a two-dimensional image (two-dimensional 

measurement data) FT, a generated three-dimensional model 
ML, other data and so on. 

The standard model DS includes a model for skin DSA, 
a model for skeleton DSB and a model for muscle 

20 (information of muscle) DSC (see Fig. 14). Common control 
points are set for the model for skin DSA, the model for 
skeleton DSB and the model for muscle DSC for a single 
standard model. Therefore, if the control points move to 
modify the model for skin DSA, for example, based on the 

25 measurement data, the model for skeleton DSB and the model 
for muscle DSC both are modified according to the movement 
of the control points. 

These programs and data are loaded on the RAM of the 
processor 10 as required. 

30 The modeling program PR includes a measurement 



process, rough alignment, a data reduction process, a 
modification process, a partial area selection process, a 
reliability obtaining process, a weighting process and 
programs for other processes. 
5 In the modification process, fitting of a surface S 

of the standard model DS is performed in order to minimize 
distances between the surface S and respective points of 
the three-dimensional data DT. Degree of the fitting 
differs according to a weight of each of points. In the 

10 weighting process, the process required for the weighting 
in the modification process such as a process of 
determining the weights using the reliability data DR. 

The media drive unit 12 is adapted to access a 
recording medium such as a DC-ROM (CD), a floppy disk FD 

15 or an optical magnetic disk to read and write data or 
program. Proper driving unit is used depending on the 
type of the recording medium. The modeling program PR can 
be installed from the recording medium. It is also 
possible to input the standard model DS, three-dimensional 

20 data DT, reliability data DR, two-dimensional image FT and 
the like via the recording medium. 

Displayed on a display screen HG of the display 13 
are above-mentioned various data, the three-dimensional 
model ML generated by the modeling program PR, other data 

25 and images. 

The key board 14 and the mouse 15 are used for 
inputting data or issuing command to the processor 10. 

The three-dimensional measuring unit 16 serves to 
obtain three-dimensional data DT of an object by, for 

30 example, the light-section method. It is possible to 



obtain three-dimensional data DT directly by the three- 
dimensional measuring unit 16, while it is also possible 
to obtain three-dimensional data DT indirectly by 
performing operations in the processor 10 or the like 
5 based on data output from the three-dimensional measuring 
unit 16 . 

In the course of obtaining three-dimensional data DT 
of an object, it is possible to obtain two-dimensional 
data FT of the object on the same visual line as that of 

10 the three-dimensional data DT when so required. A known 
device such as that disclosed in Japanese Unexamined 
Patent Publication No. 10-206132 can be used as the three- 
dimensional measuring unit 16. 

Another known method for obtaining three-dimensional 

15 data DT of an object is a method wherein a plurality of 
cameras is located with parallax for an object. 
Specifically, the three-dimensional data DT is obtained by 
computing a plurality of images with parallax captured by 
the cameras and using the stereophotography method. 

20 In the above method, data for judging reliability of 

respective points of the three-dimensional data DT can be 
obtained by using three cameras . 

More specifically, an object is captured by three 
cameras to obtain three images, and points corresponding 

25 to the three images are determined. The three-dimensional 
data DT is then obtained from a known calculation based on 
the corresponding point of two images . Another image is 
utilized for obtaining the reliability data DR. 

For example, as shown in Fig. 10, three -cameras A, B, 

30 and C capture the object Q to obtain three images FA, FB 
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and FC. Image planes (Un, Vn; n= 1, 2, 3) are indicated 
for the images FA, FB and FC. The point QP on the object 
Q in the three-dimensional space M is projected on each of 
the points PA, PB and PC of the image planes. 
5 Here, if the relationship of the points PA, PB and 

PC is determined, the three-dimensional space M' can "be 
reconstructed from the relationship. In the three- 
dimensional space M' , the point QP' which corresponds to 
the points PA and PB is determined. Logically, the point 
10 PC obtained by reversal projection of the QP' in the 

three-dimensional space M' on the image plane (U3, V3 ) is 
identical to the point PC obtained by projecting the point 
QP in the three-dimensional space M on the image plane (U3, 
V3) . 

15 However, since it is difficult to precisely 

determine the projection conversion and the relationship, 
the points typically are not identical to each other. 
Therefore, the difference between the point PC and the 
point PC is used as an index for the reliability. 

20 For example, the difference between the point PC and 

the point PC may be indicated by number of different 
pixels. If the point PC and the point PC are on the same 
pixel, the difference is "0". If the points are different 
to each other with respect to one pixel, the difference is 

2 5 "1"; and if the points are different with respect to two 
pixels, the difference is "2". The number of different 
pixels can be used as the reliability data DR. 

The reliability data DR may be judged by other 
methods such as that disclosed in Japanese Unexamined 

30 Patent Publication No. 61-125686. 
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The modeling device 1 may be provided by using a 
personal computer, a workstation or the like. It is 
possible to obtain the program and data mentioned above by 
receiving them via a network NW. 
5 Hereafter, general flow of process in the modeling 

device will be illustrated with reference to the flowchart. 
(Preparation for Standard Model) 

Referring to Fig. 2, a standard model of an object 
is prepared as a start (#11). Since the object is a human 

10 head in the present embodiment, a standard model that is 
most resembles to the object is chosen among a plurality 
of standard models having various sizes and forms to be a 
standard model DS1 . 

The standard model may be either a three-dimensional 

15 form model or a three-dimensional shape model defined by 
polygons or a three-dimensional form model or a three- 
dimensional shape model defined by free-form surfaces. In 
the case where the standard model is the three-dimensional 
form model defined by polygons, the standard model has a 

20 form depending on three-dimensional coordinates of 

vertexes of the polygons . In the case where the standard 
model is the three-dimensional form model defined by free- 
form surfaces, the standard model has a form depending on 
functions that define the free-form surfaces and 

25 coordinates of control points. 

Further, in the case where the standard model is the 
three-dimensional form model defined by polygons, the 
vertex of polygon is referred to as "construction point" 
in this specification. 

30 A point used for modifying a standard model in 



fitting the standard model is referred to as "control 
point" in this specification. Allocations of the control 
point and the construction point are arbitrary, and the 
control point may be set either on a surface of a polygon 
5 or away from a surface of a polygon. Single control point 
is associated with a plurality (about 3 to 100) of 
construction points, so that the associated construction 
points move in accordance with movement of the single 
control point. A plurality of the control points are 
10 moved to modify whole parts of the standard model in 

fitting the standard model. In the case where the three- 
dimensional form model is defined by free-form surfaces, 
allocation of control points to be used for fitting is 
arbitrary . 

15 The control points are located at a high density at 

a part having a delicate form such as a corner of an eye 
and an edge of a lip or a part having a steep and jagged 
form such as the nose and lips. The control points are 
located uniformly at other parts . 

20 On the standard model, there are set outlines and a 

characteristic point as viewed from one direction. The 
outlines are set for the eyes, nose, mouth, jaw and the 
like, respectively. Examples of the outlines include an 
eyelid outline, a nose outline, a lip outline, a jaw 

25 outline and the like. Examples of the characteristic 

points include corners of the eyes or mouth, top of the 
nose, bottom of the jaw and like parts that actually are 
characteristic and parts which are not characteristic but 
can easily be specified as being at an intermediate part 

30 of the characteristic parts. 
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Ref erring to the standard model DS1 shown in Fig. 4, 
RK1 denotes the jaw outline, RK2 denotes the lip outline 
and RK3 denotes the eyelid outline. As shown in Fig. 4, 
the outline RK1 is a boundary when the standard model DS1 
is viewed from one direction. Only a part of the 
characteristic points TT of the standard model DS1 is 
shown in Fig. 4. 

(Obtaining of Three- Dimensional Data) 

The three-dimensional data DT is obtained by 
performing three-dimensional measurements of an object 
(#12). As mentioned above, at the same time with 
obtaining the three-dimensional data, two-dimensional 
image, reliability data for indicating reliabilities of 
points of the three-dimensional data and information for 
obtaining the reliability data. 

Referring to Fig. 5, head of a person is measured 
(captured) as an object by using the three-dimensional 
measuring unit 16, to thereby obtain the three-dimensional 
data and the two-dimensional image. 

The reliability data can be obtained by the 
reliability obtaining process of the modeling program. 
The three-dimensional data and/or the two-dimensional 
image are sometimes referred to as "measurement data" in 
this specification. Order of performing the preparation 
for standard model and the obtaining of the three- 
dimensional data is not crucial, and the preparation for 
standard model and the obtaining of the three-dimensional 
data may be performed in parallel. 
(Rough Alignment) 

The standard model and the three-dimensional data 



are roughly aligned in this step (#13). In order to align 
the standard model and the three-dimensional data, 
orientation, size and location of the standard model are 
changed. The standard model is magnified in the 
directions X, Y and Z respectively at an arbitrary 
magnification, so that the size of the standard model in 
the directions approximates to that of the three- 
dimensional data. 

It is possible to obtain the standard model DSa 
shown in Fig. 6, which has nearly the same size as the 
three-dimensional data DT, by rotating and magnifying the 
standard model DS shown in Fig. 6. In Fig. 6, the 
standard model and the three-dimensional data are not 
aligned for the convenience of understanding. 

There are practically two methods of alignment, one 
is (1) general rough alignment and the other is (2) 
topical rough alignment as explained below. 

Method (1) can be performed automatically. A part 
of method (2) should be manually performed, since it is 
difficult to perform extraction of characteristic points, 
which is one of the steps included in method (2). Method 
(1) is suitable in the case where the sameness of the 
forms weighs heavily, and method (2) is suitable in the 
case where sameness of locations is rated above the 
sameness of forms such as in animation production, from 
the viewpoint that the standard model is modified 
topically in the fitting after the rough alignment. Time 
for the alignment process is effectively shortened by 
employing method (2) if a user is familiar with the 
extraction of characteristic points . 



(General Rough alignment) 

In the general rough alignment, location, 
orientation and size of the standard model is changed so 
as to minimize distance between the three-dimensional data 
5 and the standard model; i.e., there are determined si, ai 
and ti when energy function e(si, ai, ti) indicated in the 
following expression is the smallest. 

In the expression, f(si, ai, ti) is an energy 
function defined in association with the distance between 
10 the three-dimensional data and the standard model, and 
g(si) is a stabilizing energy function for avoiding an 
excessive modification. 

Initial values of the location, direction and size 
of the standard model may be obtained by performing 
15 pattern matching on two-dimensional image obtained at the 
same time with obtaining the three-dimensional data by 
means of the three-dimensional measuring unit 16. 
e(si,cci, ti)=f ( si, ai, ti) +g( si) (1) 

f(si,ai,ti)= 2 d k (si,ai,ti) 2 

k 

20 g(si)=^{(si-s 0 ) 2 +(s 2 -so) 2 +(S3-so) 2 } 

so 

wherein, 

K: number of construction points of three- 
dimensional data; 

dk: distance between construction point and standard 

2 5 model; 

Wsc: weight parameter for stabilizing magnification; 
SO: initial scale; 

Si: amount of magnification in any one of directions 
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(S3 is a depth); 

cci: rotation of standard model in any one of 
directions ; and 

ti: amount of movement of standard model in any one 
of directions . 

Here, the construction points on the standard model 
move according to the following expression (2) and 
distances between the construction points and the surface 
of the standard model change according to the movements. 

M k {sx,ai, ti,x) = Qr3Qr2QrlQsx + t (2) 

%: coordinate of point of object for modification 
t: movement vector, t=[ti,t2,t3] 
Qri: rotation matrix of %± axis. 



0 0 
cosaj sinaj 
-sin ax cosa^ 



and so on. 
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D: conversion matrix from global coordinate. 



When the three-dimensional measuring unit 16 
performs measurements (capturing) from only one direction, 

20 magnification in the depth direction (direction Z) 

sometimes cannot be precisely measured. In this case, the 
magnification in the direction Z can be corrected by using 
magnifications in directions X and Y (SI and S2), as 
indicated by the expression (3), with considering that 

25 there is not a big difference between the form of the 

three-dimensional data and the form of the standard model. 
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e(si,ai,ti)=f ( si,ai, ti)+g(si) (3) 

f (si,cci,ti) = 2 d k 2 

k 

g(si)=^{(si-so) 2 +(s2-s 0 ) 2 +Y[s 3 -^- (si+s 2 ) ] 2 } 

SQ 1 

wherein , 

5 y: weight parameter with respect to modification in 

the visual line direction x3 . 
(Topical Rough Alignment) 

In the case where the automatic alignment performed 
in the general rough alignment as mentioned above is not 

10 successful, the alignment should be performed manually. 
The topical rough alignment herein explained is a method 
wherein the manual alignment is simplified as possible. 
Unsuccessful alignment should be reset for performing the 
alignment again from the start . 

15 In the topical rough alignment, correspondences 

between the characteristic lines or points on the three- 
dimensional data and the characteristic lines or points on 
the standard model are detected, and location, direction 
and size of the standard model are changed in order to 

20 minimize distances between them. In the case where 
correspondence between two lines are detected, 
perpendicular lines are drawn from a point on one line to 
another line and points that forms the shortest 
perpendicular line are set as characteristic points. By 

25 repeating such process, a plurality of characteristic 
points is obtained on the lines . 

Specifically, ti, cci and si of the standard model is 



determined in such a manner that the energy function E(si, 
cci, ti) represented by the expression (4) becomes the 
smallest with respect to a distance between a 
characteristic point on the three-dimensional data and the 
corresponding characteristic point on the standard model. 
E(si,cci,ti)= 2 |M k (si,ai,ti,x)-C k | 2 .... (4) 

k 

wherein , 

k: number of corresponding characteristic points; 
Mfc: characteristic points on standard model after 
alignment ; 

x: characteristic points on standard model before 
alignment ; 

C k : characteristic points on three-dimensional data; 
si: magnification of standard model in each 
direction; 

ai: rotation of standard model in each direction; 

and 

ti: movement of standard model in each direction. 

When the three-dimensional measuring unit 16 
performs measurements (capturing) from only one direction, 
magnification in the depth direction (direction Z) 
sometimes cannot be precisely measured. In this case, the 
scale in the direction Z can be corrected in the same 
manner as in the general rough alignment, as indicated by 
the following expression (5). 



E(si,cti,ti) = 

2 |M k (s±,ai,ti,x)-C k | 2 + Yts 3 -|(si + s 2 }3 2 (5) 

k 2 

(Extraction of Outlines and Characteristic Point) 

In this step, outlines and characteristic points are 
extracted on three-dimensional data DT or a two- 
dimensional image FT (#14). The outlines and 
characteristic points of a standard model can be extracted 
in advance of the step #14. In such cases, the outlines 
and the characteristic points which should be located at 
the same locations as those of the outlines and 
characteristic points on the standard model are located on 
three-dimensional data DT or a two-dimensional image 
corresponding to the three-dimensional data DT (see Fig. 
7). 

In the case where the outlines and characteristic 
points are not extracted in advance, they are assigned on 
the standard model at the same time when they are located 
on the three-dimensional data or the two-dimensional image. 
(Data Reduction) 

Next, in order to reduce number of calculations and 
errors, reduction of the three-dimensional data is 
performed by extracting data that are necessary and highly 
reliable (#15). This data reduction enables to reduce the 
number of calculations without deforming the form of the 
original three-dimensional data. 

Unnecessary data such as data for outside an object 
are eliminated in the data reduction. For example, an 
area for the face of a two-dimensional image is judged in 



order to maintain only the three-dimensional data 
corresponding to the area. The area can be judged using 
differences in the distances from background of the object 
and parts of the object. The area of the face can be 
judged by using the rough alignment information. In the 
case where the reliability data accompanies the three- 
dimensional data, only the data of high reliability are 
maintained. When data density is high, some of the data 
are deleted so that the density is uniform. 

In the case of deleting data for making the density 
uniform, only the data satisfying the following expression 
(6) are maintained. 

|Pk-P-r|>^ et(Pk)+ 2 ^et( p ^) (6) 

wherein , 

Pk: construction points; 

P-r: construction points used previously; and 
Rdet (x) : function which represents density of 

periphery of a construction point x. 

According to the above expression (6), the data Pk 

is used if the distance between the data Pr used 

previously and the data Pk is more than a predetermined 

value . 

(Modification) 

In this step, modification of whole of the standard 
model DS is performed (#16). There are used, to determine 
a general energy function e, an energy function el which 
is defined in association with the distances between the 
construction points of the three-dimensional data and the 
surface of the standard model; an energy function e3 which 
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is defined in association with the distances between 
characteristic points of the standard model and the 
characteristic points assigned for the three-dimensional 
data; an energy function e2 which is defined in 
association with the distances between the outlines of the 
standard model and the outlines assigned for the three- 
dimensional data; and an energy function es which is 
defined to avoid an excessive modification. The general 
energy function e determined is evaluated. The surface of 
the standard model is modified in such a manner that the 
general energy function e becomes the smallest value. 

It is preferred to use el, e2, e3 and es as a 
general energy function e; however, it is possible to use 
arbitrarily two functions among el-e3. 

Each of the energy functions will be explained in 
the following. 

(Distance Between Standard Model and Three -Dimensional 
Data) 

Referring to Fig. 8, the point Pk denotes a point in 
a point group that is comprised in the three-dimensional 
data DT. The point Qk indicated in the surface S of the 
standard model DS is the nearest point from the point Pk. 
The point Qk is an intersection point that is found by 
drawing perpendiculars from the point Pk to the surface S. 
The distance between the point Pk and the point Qk is 
evaluated in this step. 

The difference energy el, which represents 
differences between the points in the three-dimensional 
data and the surface of the standard model, is determined 
by the following expression (7) using a square distance of 



the distance between the point Pk on the three-dimensional 
data after the data reduction and the point Qk found by 
projecting the point Pk on the surface S of the standard 
model DS. 

el(TlA) = -^2^)[d kT (Qk-Pk)] 2 (7) 

wherein, 

T1A: control points; 

Pk: construction point of three-dimensional data 
after data reduction; 

Qk: projecting point of the construction point of 
three-dimensional data on the surface of model; 

K: number of construction points after data 
reduction; 

dk: direction of projection from construction point 
to surface of model: d k = (Qk-Pk) /|Qk-Pk|; 

pk: reliability of construction point Pk; 
w(pk): reliability function, w(pk)=l/(a + pk) n ; 
W: 2w(pk); and 

L: adjusting scale that enables to handle various 
energy by one common unit . 

{Distance between Outline of Standard Model and Outline of 
Measured Data) 

In this step, there are evaluated distance between 
an outline assigned on three-dimensional data or an 
outline assigned on two-dimensional image and an outline 
on a standard model. 

When the outline of the measured data is assigned on 
the three-dimensional data, the distance therebetween is 
determined by drawing perpendiculars from points on the 
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outline of the three-dimensional model to the outline on 
the standard model which corresponds to the outline of the 
three-dimensional data and selecting the shortest 
perpendicular as the distance. 

In the case where the outline of the measured data 
is assigned on the two-dimensional image, the outline of 
the standard model is projected on the two-dimensional 
image by using camera parameters for the camera which has 
captured the two-dimensional image. Distance is 
determined by drawing perpendiculars from points on the 
outline of the two-dimensional image to the outline of the 
standard model that corresponds to the outline of the two- 
dimensional image and selecting the shortest perpendicular 
as the distance. 

In the case where the outline of the measured data 
is assigned on the three-dimensional data, e2 which 
represents difference energy of respective outlines of the 
standard model is calculated by the following expression 
(8) using sum of squares of the distances. 

e2(T2A) = -Vi [d kT (gk-pk) ] 2 (8) 

nl 2 kti 

wherein, 

T2A: control points; 

pk: points on outline on three-dimensional data; 

qk: pedal point of perpendicular line drawn from 
point on outline of three-dimensional data to 
corresponding model outline; 

n: number of points on outline of three-dimensional 
data associated with one model outline; 

d : direction of projection of outline of measured 
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data to corresponding model outline, 
d k = (qk-pk)/|gk-pk|; and 

1: adjusting scale for handling various types of 
energy by one common unit . 

In the case where the outline of the measured data 
is assigned on the two-dimensional image, e2 ' which is the 
difference energy of respective outlines of the standard 
model is calculated by the following expression (9). 

e2'(T2'A) = J- | [d kT (gk-pk)] 2 (9) 

nl 2 jfeti 

wherein , 

T2A: control points; 

pk: points on outline of two-dimensional image; 

qk: pedal point of perpendicular line drawn from 
point on outline of two-dimensional image to corresponding 
model outline projected on two-dimensional image; 

n: number of points of outline of measured data 
associated with one model outline; 

dk: direction of projection from points on outline 
of two-dimensional image to corresponding model outline, 

d k = (qk-pk)/|qk-pk| ; and 

1: adjusting scale for handling various types of 
energy by one common unit . 

Reason for assigning the outline of the measured 
data on the two-dimensional image is that the outline 
itself can be inaccurate if the outline is assigned on the 
three-dimensional data that are ambiguous. Thus, the 
outline is extracted using the two-dimensional image. 
(Distance between Characteristic Points on Standard Model 
and Corresponding Characteristic Points on Measured Data) . 



Characteristic points are set on the measured data 
in order to evaluate distances between the characteristic 
points on the thee -dimensional data or the characteristic 
points assigned on the two-dimensional image and the 
characteristic points of the standard model. 

Difference energy e3 between the characteristic 
points on the three-dimensional data and the 
characteristic points on the standard model is calculated 
from the following expression (10) by using a square 
distance of the distance between corresponding 
-characteristic points. 

In the case where the characteristic points are 
assigned on the two-dimensional image, the difference 
energy on the two-dimensional image is calculated by 
projecting the characteristic points of the standard model 
on the two dimensional image using camera parameters . 

e3(T3A) = -i-f |F k -G k }| 2 .... (10) 
NlS jfc.i 

wherein , 

T3A: control points; 

Fk: characteristic points of measured data; 

Gk: characteristic points on standard model 
corresponding to characteristic points of measured data; 

N: number of correspondences among characteristic 
points on standard model and characteristic points of 
measured data; and 

L: adjusting scale for handling various types of 
energy by one common unit. 

(Stabilization Energy for Avoiding Excessive Modification) 
In addition to the difference energy mentioned above. 



a stabilization energy es is used for avoiding an 
excessive modification. 

Specifically, distances between control points used 
for modification are connected by virtual springs (elastic 
bars) KB as shown in Fig. 9. Based on restrictions of the 
virtual springs KB, the stabilization energy es for 
stabilizing the form of the surface S of the standard 
model DS is defined. 

The virtual springs do not necessarily connect the 
control points, but relationship between the control 
points and the virtual springs must be definite. 

A part of the surface S of the standard model DS, 
which is an object for fitting, is shown in Fig. 9. The 
surface S is formed of control points U = |ui, i = 1 . ...n|. 
A virtual spring KB is provided between neighboring 
control points . The virtual spring KB restrains the 
control points by tensile force and serves to prevent the 
surface S from improper modification. 

Specifically, the tensile force of the virtual 
spring KB is increased with the increase in an interval 
between neighboring control points U. For example, the 
tensile force of the virtual spring KB increases when the 
point Qk approaches the point Pk, and the interval between 
the control points U increases according to the movement 
of the point Qk. The tensile force of the virtual spring 
KB does not change when the movement of the point Qk does 
hot affect the interval of the control points U, i.e., 
when relative position relationship between the control 
points U does not change. The stabilization energy es is 
defined as an average tensile force of the virtual springs 



-28- 
KB on the whole surface S. Accordingly, the stabilization 
es increases when a part of the surface S is excessively- 
modified. The stabilization es is null when the surface S 
is uniformly modified. 
5 The stabilization energy es is determined by the 

following expression (11) according to the state of 
deformation of the virtual spring KB. 

es(TsA) = -^j 2 t-^j- (U~m-V~m) T (Um-Vm)-L0m] 2 (11) 

ML m=l Lq" 

wherein , 
10 T3A : control points; 

U~m, V~m: initial value of edge (control point) of 
virtual spring; 

Urn, Vm: edge of virtual spring after modification; 
LOm: length of virtual spring in initial state, 
15 LOm = |U~m - V~m( 

M: number of virtual springs; 
c: spring coefficient; and 

L: adjusting scale for handling various types of 
energy by one common unit. 
20 Therefore, the virtual spring KB is hardly deformed 

if the spring coefficient c is large. 

The stabilization energy function es serves as a 
certain restraint for modification of the surface S and 
thus enables to prevent the surface S from an excessive 
25 modification. 

(General Energy Function) 

As described in the foregoing, the control points 
T1A, T2A, T3A and TsA are used for the energy functions el, 
e2, e3 and e4 , respectively. Although the identical 



control points TIA-TsA are used in the foregoing, control 
points different from one another can be used as mentioned 
below. The standard model is modified using the different 
control points TA, and the control points TA as 
represented by the expression (12) which can minimize the 
general energy function e(TA) is determined. 
e(TA) = 

wlel(TA) + 2 w2e2 S (TA)-fw3e3(TA)+ces(TA) (12) 

wherein , 

el(TlA): difference energy between construction 
points of three-dimensional data and surface of model; 

e2 s (T2A): difference energy between model outline 
and outline of measured data; 

e3(T3A): difference energy between characteristic 
point of measured dta and characteristic point of model; 

eS(TSA) : stabilization energy for avoiding excessive 
modification; 

wi, c: weight parameter for respective energy; and 
TA = T1A = T2A = T3A = TSA. 
(Repetitive Modification) 

In actuality, the modification is performed 
repeatedly (#17). Specifically, the same step as the step 
16 is repeated after changing control points to perform 
modification. The modification can be performed by 
changing the energy functions instead of changing the 
control points. Assuming that the general energy function 
at n time(s) of modification is en(TA), the general energy 
function en(TA) is considered to be converged when the 
following expression is satisfied. 
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|en(TA) - en-l(TA)| < e (13) 

Hereinafter, a general modification process is 
explained with reference to Fig. 3. Firstly, a pair of 
5 corresponding points between the measured data and the 
standard model is prepared (#21) . For example, the point 
Qk that corresponds to the point Pk shown in Fig. 8 is 
calculated and determined to prepare a pair of the points 
Pk and Qk. 

10 The surface S is modified (#22), and the general 

energy function en(TA) after the modification is 
calculated (#23). Until the general energy function 
en(TA) is converged (until it is "YES" in #24), the 
process is repeated. 

15 There may be employed a known method for judging the 

convergence of the general energy function en(TA) such as 
a method wherein the convergence is confirmed when the 
general energy function en(TA) is smaller than a 
predetermined value and a method wherein the convergence 

20 is confirmed when degree of change from the previous 
calculation is smaller than a predetermined value. 
(Different Control Points) 

In the expression (12), the same control points are 
used for the different energy functions el-e4 having 

25 different objects for fitting (e.g., construction points 
of three-dimensional data, outline, characteristic point). 
In turn, control points are varied for each object for 
fitting, i.e., for each energy functions in the following 
example. Specifically, in the following example, 

30 different control points T1A, T2A, T3A and TSA are used. 



In this case, the expression (12) can be used as the 
general energy function e(TA) as long as the control 
points T1A, T2A, T3A and TSA are different from one 
another and in the following relationship. 
5 TA D T1A 

TA D T2A 

TA D T3A 

TA = TSA 

As mentioned above, movement due to energy between 
10 the control points is topical in nature since the control 
points are topical (as points usually are) . For example, 
in aligning the position of an eye of three-dimensional 
data and the position of an eye of a standard model, 
undesired modification may occur if a part where 
15 characteristic points are set is strongly pulled. Such 
undesirable modification should be avoided. 

In turn, the construction points of the three- 
dimensional data need to be located to realize delicate 
movement. However, in the case where the control points 
20 are scarcely set, the construction points cannot move 
delicately. 

Therefore, it is preferred to use a large number of 
control points with respect to the construction points of 
the three-dimensional data and a small number of control 
25 points with respect to the characteristic points. An 
intermediate number of control points are used for the 
outline. 

The control points are set with high density for an 
outline that has a steep or jagged form. The 
30 stabilization energy must be applied to all the control 
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points. The selection of control points as mentioned 
above is performed in preparation for the standard model. 

The control points T1A, T2A, T3A and TSA are 
different from one another; however, control points 
included in each of the T1A, T2A, T3A and TSA can be used 
in common. 

(Change in Weight Depending on Reliability) 

In the expression (12), the general energy function 
e(TA) is evaluated assuming that reliabilities of 
information are equal to one another; however, in this 
example, weights are changed depending on the 
reliabilities of the information. Thus, it is possible to 
judge e(TA) with counting strongly the information of high 
reliability. 

The reliability of respective information can be 
obtained in the three-dimensional measurement or automatic 
extraction of outlines and characteristic points. 

In this case, control points (TA) which can minimize 
the general energy function e(TA) represented by the 
expression (14) is determined. 

e(TA) = wlWKpl)el(TA) + £ w2W2(p2s)e2s(TA) 

+ w3W3(p3)e3(TA)+ces(TA) (14) 

wherein , 

pi: reliabilities of information regarding 
respective energy functions ei(TA); 

W(i): reliability function; and 
TA = T1A = T2A = T3A = TSA. 

Fitting with weighting based on the reliability will 
hereinafter be described in more detail. 



-33- 

Fitting energy function Ef (U) for a point Pk in a 
point group, a point Qk corresponding to the point Pk and 
corresponding points T = {(Pk, Qk) , K = 1 . . . .n} may be 
set as in the following expression (15), for example. 

Ef(U) = | ||Pk-Qk(U)|| 2 (15) 

wherein, Qk(U) indicates that Qk is a function of U. 

Stabilization energy function Es(U) is represented 
by the following expression (16). The expression (16) 
corresponds to the expression (11). 

c M •( m 

Es < u ) =^2 t — (n~r-n~-) (nf-n-)-L-] .... (16) 

m-l Lq 

Evaluation function E(U) of fitting is defined as 
the following expression (17) using the fitting energy 
function Ef(U) and the stabilization energy function Es(U). 

E(U) = WfEf(U) + WsEs(U) (17) 

In the above expression, each of Wf and Ws is a 
weight coefficient used for normalization. 

The fitting of a surface is performed by repeating 
modification of the surface and search for corresponding 
point so that the evaluation function E(U) of the 
expression (17) becomes satisfactorily small. For example, 
the surface is moved in the fitting to the direction in 
which differential of U of the evaluation function E(U) 
approaches to null. 

The evaluation function E(U) described above is 
obtained by assuming that the reliabilities of the points 
are equal to one another and handling the reliabilities of 
points of a point group as equal to one another. 
Weighting based on the reliability of the point Pk is 



conducted with respect to the evaluation function E(U) . 

Specifically, in the case where the fitting energy 
function represented by the expression (15) is used, an 
undesirable fitting result may sometimes be introduced as 
5 being affected by the points of low reliability since the 
data for all the points Pk cannot be fully reliable. 
However, in the case where the point Pk of low reliability 
is deleted, a point group in which the point Pk is deleted 
may not be able to maintain the original form thereof. 
10 Therefore, elements of the reliabilities of points 

Pk in a point group is introduced to the expression (15) 
as weights of corresponding relationships to obtain the 
following expression (15'). 

Ef(U) = 2 w(pk)||Pk-Qk(U)|| 2 (15') 

15 In the above expression (15'), pk represents the 

reliability of the point Pk. It is possible to use the 
above-mentioned reliability data as the reliability pk. 
W(pk) is a function for determining the weight from the 
reliability pk. 

20 Generally, in the case where x represents 

reliability of a point, W(x) is a function which serves to 
assign a weight to the point depending on its reliability. 
For example, large weight is assigned if the reliability x 
is high (when value of x is small), and small weight is 

25 assigned if the reliability x is low (when value of x is 
large) . Such function can be obtained by the following 
expression (18). 

W(x) = l/(l+x) n (18) 

In the expression (18), n is a weighting coefficient 



for the reliability, and x > 0. 

For example, in the case of using number of pixels 
which is the reliability data as the reliability x, x = 0 
and the weight W(x) = 1 if there is no shift. When x = 1, 
5 the weight W(x) is l/2 n ; and when x = 2, the weight W(x) 
is l/3 n . When the reliability decreases (when value of x 
increases), the weight W(x) approaches to null. Thus, the 
surface is modified in accordance with the reliability of 
each of the points Pk. 

10 When the expression (15') is used as a fitting 

energy function Ef(U) in the expression (17), more precise 
fitting is performed with respect to the point Pk of high 
reliability and fitting with respect to the point of low 
reliability is performed is such a manner as to avoid 

15 affecting to the fitting with respect to the point of high 
reliability. 

For example, referring to Figs. 18A and 18B, in the 
case where the points PI and P2 have high reliabilities 
and points 3 to 6 have low reliabilities among the six 

20 points of points 1 to 6, the surface S is modified in 

accordance with the reliabilities of the points 1 to 6 . 
In the surfaces Sc and Sd modified according to the 
present embodiment, the form formed of a point group 
reflects the reliabilities of the points, the form of the 

25 point group is maintained and excessive modification is 

suppressed as compared to the surface Sa in Fig. 19A which 
is modified without considering the reliabilities of the 
points and the surface Sb in Fig. 19B which is modified 
with deleting the points of low reliability. 

30 As described above, the degree of fitting of a 



surface to the point group comprising a set of points Pk 
is changed based on the reliability information of the 
points Pk according to the present embodiment. The 
present embodiment thus prevents improper modification or 
erroneous modification even when points of low 
reliabilities are included in a point group, thereby to 
achieve a proper fitting. 

Further, since the modification of surface is 
performed depending on reliabilities of points Pk without 
deleting points of low reliabilities as mentioned above, 
the modification reflects an original form of a point 
group. Thus, the present invention achieves more proper 
fitting than that of the conventional method of deleting 
the points of low reliabilities. 

Although reliabilities of points of a point group is 
judged and then a weight of each of the points is changed 
based on the Judgment in the present embodiment, it is 
also possible to judge the reliabilities by a 
predetermined area comprised in a point group and then 
change weights of all the points existing in the area at 
once based on the judgment . 
{Partial Area Modification) 

By the steps #1 to #17, it is possible to obtain a 
three-dimensional model more elaborated and precise than 
that obtained by the conventional method. However, 
process of the steps #1 to #17 is performed with respect 
to the whole part of an object and, therefore, a fitting 
with respect to a part of the object having topical 
characteristics may not sometimes be satisfactorily 
precise. Accordingly, the characteristic part of the 
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object or a part with respect to which a fitting have not 
been performed satisfactorily in the steps #1 to #17 is 
extracted from three-dimensional data, and a fitting is 
performed using the extracted three-dimensional data (#18). 
5 Fig. 11 is an example of the extraction of partial 

areas BRY, and Fig. 12 is an example of a method for 
modifying a standard model DS using the partial areas BRY. 

Referring to Fig. 11, the partial areas BRY are 
extracted from three-dimensional data DT . The partial 

10 areas BRYa, BRYb and BRYc are three-dimensional data of 
the eyes and the periphery, the mouth and the periphery 
and the nose and the periphery, respectively, which are 
extracted from three-dimensional data DT of an object. 
The partial areas BRY may be reduced in data when so 

15 required. The data reduction can be performed in the same 
manner as that of the step #15. Since fitting area is 
limited to the partial areas BRY in a modification of 
partial area described below, the right-hand side of the 
expression (6) is made closer to null than that of the 

20 step #15 to lower the data reduction rate; therefore, the 
fitting process does not take time even in the case of 
maintaining a certain level of precision of the partial 
areas BRY. Therefore, it is possible to set the data 
reduction rate of each of the partial areas BRY 

25 considering a region of the partial area BRY, required 

precision and process time. The extraction of the partial 
areas BRY may be conducted either in the step #12, i.e., 
when obtaining the three-dimensional data or after the 
process of the steps #1 to #17 is completed. 

30 Referring to Fig. 12, process of partial area 
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modif ication of the standard model DS1 indicated at the 
left-hand side is shown in the standard models DSla and 
DSlb at the right-hand side. 

In the step #18, the standard model DS1 obtained by 
5 the steps #1 to #17 is subjected to fitting using the BRYa 
as shown in the center of Fig. 12 in order to obtain the 
standard model DSla. The fitting in the step #18 is 
performed in the same manner as that of the step #16, i.e., 
by modifying the standard model and repeating the 

10 modification. In the step #18, there may be settled, 

depending on the partial areas BRY, a suitable value for 
each of the control points T1A, the point Pk and the Point 
Qk in the expression (7); the control points T2A, the 
point pk and the point qk in the expressions ( 8 ) and ( 9 ) ; 

15 the control points T3A, the point Fk and the Point Gk in 

the expression (10); the control points TsA, the point U~m 
and the Point V~m in the expression (11); the threshold 
value £ of the expression (13) or the like. 

The partial area modification is performed 

20 successively for each of the partial areas BRY (#18 and 
"YES" in the step #19). Control points and the like are 
varied for each of the partial areas BRY, so that fitting 
of the standard model DS1 can be more precisely performed. 
For example, after obtaining the standard model DSla, 

25 fitting is performed using the partial area BRYb to obtain 
the standard model DSlb. The standard model DSlb thus 
obtained is the target three-dimensional model ML. 

Each of the partial areas BRY may be overlapped with 
a part of other partial areas BRY or included in any of 

30 other partial areas BRY. For example, referring to Fig. 



-39- 

11, a part of the partial area BRYa may be overlapped with 
a part of the partial area BRYc , or a partial area 
modification of the partial area BRYb may be performed 
followed by an extraction of data of the upper lip as a 
new partial area and then a partial area modification for 
the new partial area. 

The extraction of partial areas may be conducted 
either by an operator (manually) or by performing the 
partial area selection process in the modeling program PR 
as shown in Fig. l. In the latter instance, an area 
corresponding to a part to be extracted as a partial area 
is defined in a standard model in advance of the 
extraction. Then, an area of three-dimensional data 
corresponding to the standard model modified by the 
process of the steps # 11 to #17 is judged as the partial 
area to be extracted. 

According to the embodiment described above, the 
fitting of a standard model is followed by the partial 
fitting of the standard model using the partial areas, so 
that fitting of topical areas such as the eyes or the 
mouth can be performed more precisely without causing 
improper modification. Thus, it is possible to generate a 
three-dimensional model that closely resembles to an 
object . 

(Generation of Skeleton Model) 

The model for skin of a standard model is modified 
to generate a skin model of a three-dimensional model by 
the above process. Referring to Figs. 13 to 17, a method 
for generating a skeleton model MLB by modifying a model 
for skeleton DSB will be explained below. 
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As mentioned above, the model for skin, the model 
for skeleton and the model for muscle of a standard model 
have control points in common. Accordingly, the control 
points set on the standard model are changed if the model 
for skin DSA is modified for generating the skin model MLA. 
The control points of the skin model MLA, i.e., the 
control points after modification are the modification 
information or modification result in the present 
invention. 

(Generation of Skeleton Model by Adopting Modification 
Information) 

The model for skeleton DSB is modified using 
modification information to generate the skeleton model 
MLB as shown in Figs. 13 to 17. 

Specifically, the model for skeleton DSB is modified 
by using the control points ThA set on the standard model 
DS and the control points TA after modification. Here, 
"h" in the abbreviation "ThA" denotes "hat" and indicates 
that the model is not yet modified. This is applicable to 
the abbreviations appearing in the following such as "VhQ" 
and "VhR". 

Point group (construction point group) Vhg forming 
the model for skin DSA before modification can be 
represented by the following expression (19) as a linear 
joint of the control points ThA. 

F Q (T A ) = ag r A .... (19) 

wherein, VhQ: point group forming model for skin 

DSA; 

aAQ: coefficient that depends on each point of point 
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groups BhQ and ThA, 

ccAQ = fA(VhQ); and 

ThA: control points before modification. 

The construction point group VQ of the model for 
skin DSA after the modification is represented by the 
following expression (20). 

V Q (T A ) = a^T* (20) 

The construction point group VhR of the model for 
skeleton DSB before modification is represented by the 
following expression (21) . 

K R (T A ) = aJr A (21) 

wherein, VhR: point group for forming model for 
skeleton DSB; 

aAR: coefficient that depends on each of point 
groups VhA and ThA, 

aAR = fA(VhR); and 

ThA: control points before modification. 

The coefficient aAR is determined by the above 
expression (21). 

The construction point group VR of the model for 
skeleton DSB after modification (skeleton model MLB) is 
represented by the following expression (22). 

V R (T A ) = o^T A (22) 

The coefficient aAR determined above and the control 
point group TA of the model for skin DSA after 
modification are substituted with the right-hand side of 
the expression (22). This substitution enables to obtain 
the construction point group VR of the model for skeleton 



DSB after modification by that very simple operation and 
thus eliminates the complicated operation as performed in 
the modification of the model for skin DSA. 
(Correction in Adoption of Modification Information) 

The modification information (i.e., the control 
point group TA after modification) is taken out of the 
model for skin DSA after modification and used for 
modifying the model for skeleton DSB in order to generate 
the skeleton model MLB as described above. In this case, 
improper modification may occur when then modification 
information is adopted as it is since the model for skin 
DSA and the model for skeleton DSB are not identical. 

For example, referring to Fig. 13, the standard 
model DS comprises the model for skin DSA and the model 
for skeleton DSB; the model for skin DSA has many 
construction points Q; and the model for skeleton DSB has 
many construction points R. Further, the model for skin 
DSA and the model for skeleton DSB have control points U 
(A, B, C and D) in common. 

Here, only three control points U which are closest 
respectively to the construction points Q and R are used 
for each of the model for skin DSA and the model for 
skeleton DSB. Accordingly, each of the construction 
points is controlled by the three closest control points 
and not influenced by any other control points. 
Specifically, for example, each of the coefficient aAQ and 
ocAR of any other control points is set as null. 

The control points U (A, B, C) are used for the 
construction point Ql of the model for skin DSA. However, 
the construction point R2 of the model for skeleton 
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corresponding to the construction point Ql of the model 
for skin DSA is not the control points U (A, B, C) , but 
control points U (A, B, D) . Proper modification may not 
be performed if such control points are used. 

For example, in the case where the control points 
U(B) move upwards and control points U(D) move downwards, 
the construction point Ql of the model for skin moves 
upwards but the construction point R2 of the model for 
skeleton DSB moves downwards. However, the construction 
point Ql and R2 are corresponding points for the skin and 
skeleton and it is necessary that they move in the same 
direction. The irregular movement may deform, the form of 
the model. 

To avoid such deformation, correction of control 
points U is performed with respect to each of the 
construction points that are objects for modification. 
The correction may be performed by the following two 
methods . 

One of the methods employs the construction point Q 
on the model for skin DSA that is closest to the 
construction point Q of the model for skeleton instead of 
using the construction point R for determining the 
coefficient a as shown in Fig. 14. 

Specifically, in Fig. 14, the construction point Ql 
on the model for skin DSA closest to the construction 
point Ql of the model for skeleton DSB is used instead of 
using the construction point R2 . By setting a certain 
correspondence between the construction point R2 and 
construction point Ql, it is possible to use the control 
points U of the construction point Ql as the control pints 



U of the construction point R2 . In this case, the 
distance between the construction point R2 and the 
construction point Ql is ignored. 

The same control points are thus used for the model 
5 for skin DSA and the mode for skeleton DSB , so that the 
improper modification mentioned above can be prevented. 

Another method employs a point RQ, instead of using 
the construction point R of the model for skeleton DSB, 
which is obtained by projecting the construction point R 
10 on the skin model MLA. 

Specifically, the modification information is 
adopted to the model for skeleton DSB. In this case, the 
control point group VhR of the model for skeleton DSB 
before modification is represented by the following 
15 expression (23) instead of the expression (21). 

V RQ (T A ) - ccf T A .... (23) 

wherein, VhRQ: point group obtained by projecting 
point group forming model for skeleton DSB on model for 
skin MLA; 

20 aARQ: coefficient depending on each points of point 
group VhR, ThA and point group VhQ, 
aARQ = f A(VhRQ) ; and 
ThA: control point group before modification. 
The coefficient aARQ is determined by the expression 
25 (23) and a construction point group VR of the model for 
skeleton DSB after modification (skeleton model MLB) is 
determined by the following expression (24) . 
V*(T A ) = af T A (24) 

In the example shown in Fig. 15, the construction 
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point R2 of the model for skeleton DSB is projected on the 
skin model MLA. Specifically, perpendiculars is drawn 
from the construction point R2 of the model for skeleton 
DSB to the skin model MLA, and an intersection (nodal 
5 point) of the shortest perpendicular is determined as a 
virtual construction point RQ2 ' . Control points U for the 
virtual construction point RQ2 ' are used as the control 
points U for the construction point R2 . 

Such correction can prevent an unwanted modification. 

10 Following is an explanation for general process of 

generating a three-dimensional model. 

Referring to Figs. 16 and 17, standard model DS 
containing at least two models is prepared (#1). In this 
step, each of the models has control points u in common. 

15 The model for skin DSA is fitted to three- 

dimensional data and modified in order to be used as skin 
model MLA (#2). Control points u are then extracted from 
the skin model MLA (#3) . Model for skeleton DSB is 
modified by using the control points u (#4) . The control 

20 points u used for the modification of the model for 

skeleton DSB have been corrected before the modification. 
The skeleton model MLB and the skin model MLA thus 
obtained are combined to form the three-dimensional model 
ML. 

2 5 The skeleton model MLB is generated by using the 

modification information of the skin model MLA in the 
above example, and the same method is applicable to 
generation of a muscle model MLC. Further, the above- 
mentioned method may be used for modifying two models 

30 corresponding to each other. 



More specifically, in the case where the standard 
model DS comprises a first model and a second model, the 
first model may be subjected to a fitting to three- 
dimensional data so as to be modified (first step) , and 
5 modification information obtained by the modification in 
the first step is used for modifying the second model 
( second step) . 

If the standard model comprises three-dimensional 
information such as muscle information, the same 

10 modification as mentioned above can be performed. Muscle 
is formed of lines that are arranged three - dimension ally 
and, therefore, the muscle information is information 
comprising three-dimensional lines. 

According to the above embodiment, a three- 

15 dimensional model such as a skin model and a skeleton 

model is obtainable by performing a fitting using three- 
dimensional data obtained by measuring an object and 
modifying a model for skin and a model for skeleton. It 
is possible to prevent abnormal modification from 

20 occurring by correcting modification information in this 
embodiment . 

Further, in the above embodiment, the configuration 
of the modeling unit 1, the circuit, the contents of 
process, the order of process, the timing of process, the 
25 setting of coefficient and the like can be changed 
according to the scope and spirits of the present 
invention. 



